Point defects on III-V semiconductor surfaces 
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I. INTRODUCTION 

The role of intrinsic point defects in determining the 
electronic structure of surfaces is a topic that has been 
debated extensively. They may act as compensation cen- 
ters and may thus be responsible for Fermi-level pinning 
at the surface. Consequently, electrically active point de- 
fects have been commonly assumed to determine Schott- 
ky-barrier heights on metal-covered surfaces . Point 
defects might also act as nucleation centers for crys- 
tal growth and thus influence the surface and interface 
morphology of devices. Despite the importance of these 
issues little is known about the basic properties of sur- 
face point defects, such as atomic geometry, position of 
charge-transfer levels, and equilibrium concentration. In 
the following we will describe how first-principles calcula- 
tions can be employed to identify the basic properties of 
surface point defects. Calculations have been performed 
for GaAs, GaP, and InP. In the present paper we will 
limit ourselves mainly to GaAs as a model system. Re- 
sults for vacancies on GaP || and InP |1 were published 
recently. 



II. COMPUTATIONAL METHOD 

In order to investigate the above noted questions, we 
performed density-functional theory (DFT) calculations 
employing the local-density approximation (LDA) for ex- 
change and correlation. Details about the method and 
extensive convergence tests can be found in Refs. P,||. 
Surface point defects were calculated with six-layer thick 
slabs. The (110) surface had a (2 x 4) periodicity in 
the [001] and [110] direction, respectively (for the defi- 
nition of lattice orientations see the coordinate system in 
Fig. |l]a). Test calculations showed that these supercells 
are large enough to calculate relaxation geometries, elec- 
tronic structure, and formation energies of point defects 
accurately 

Using the total energy E as calculated for the defect 
and the surface systems we obtain the defect-formation 
energies Ef as a function of the cation or anion chemi- 
cal potential and the position of the surface Fermi level 



-Epermi @- Here we chose the anion chemical potential 
/i As- To give an example, the formation energy of an 
arsenic surface vacancy in the charge state q, V£ s , is: 

E f( V D = E(VD - £( 8lab ) + MAs + g^Fcrmi . (1) 

Here E(V£ S ) is the total energy of the system with and 
E(s\ab) that of the system without a vacancy. The chem- 
ical potential of As is controlled by the As partial pres- 
sure and temperature. It can be varied from gallium-rich 
(^As = —AHf) to arsenic-rich conditions (,uas = 0) [Q, 
where AHf = 0.7 eV is the heat of formation of GaAs. 
In Fig. E{ is plotted as a function of £^Fermi- For each 
point defect and i?Fcrmi only the charge state with mini- 
mum formation energy is shown. According to Eq. (|l|) 
the charge state of a defect is given by the slope of the 
curve. The points where the curve changes it's slope de- 
fine the positions of i^Fermi at which the defect changes 
it's charge state, i.e. the charge-transfer levels of the 
defect. 



III. ATOMIC AND ELECTRONIC STRUCTURE 
OF POINT DEFECTS 

For several reasons the (110) surface of III-V semicon- 
ductors is ideally suited for the study of intrinsic point 
defects, (i) It is the cleavage plane of these materials and 
high quality surfaces with a low concentration of steps 
and defects can be prepared by cleaving wafers. Scan- 
ning tunneling microscope (STM) experiments achiev- 
ing atomic resolution on these surfaces showed point de- 
fects such as vacancies and adatoms in a detailed man- 
ner [p|-pd]| . (ii) The surface does not reconstruct but un- 
dergoes a relaxation that is characterized by a buckling 
of surface anions and cations. We calculate the buckling 
angle to be 29° for GaAs in good agreement with DFT 
calculations JlJ] as well as experimental low-energy elec- 
tron diffraction analysis ||l3|| . The point-group symmetry 
of the surface is Cih with a single mirror plane pointing 
in the [110] direction, (iii) With the exception of GaP 
the fundamental band gap of the defect-free surface is 
free from surface bands. Therefore the presence of point 
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FIG. 1. Atomic relaxations of surface point defects on 
GaAs(llO): antisites (a, b), vacancies (c, d), and adatoms 
(e, f). 



defects with energy levels in the band gap might notably 
affect the electronic properties of the surface. 



A. Antisite defects 

Our calculations show that arsenic antisite Asc a de- 
fects in the surface layer of GaAs show a rather large 
displacement of 0.5 A in the [001] and [110] directions 
(Fig.|l]a). Thus the defect atom is nearly in the same 
plane as the surface anions: The large surface buckling 
between anions and cations of the defect-free surface is 
locally lifted. The nearest neighbor atoms of AsGa and 
all other atoms in the cell show only small relaxations 
below 0.2 A and 0.05 A, respectively. This rather local- 
ized relaxation pattern is a typical feature of all surface 
point defects investigated. The atomic structure of the 
defect-free surface remains nearly unperturbed beyond 
second nearest neighbor atoms of surface point defects. 
As a further common feature of most point defects we 
find that the symmetry of the surface is well conserved 
by the defect relaxation. The only exception where we 
calculate a nonsymmctric relaxation to have a preferable 
formation energy is the positively charged anion vacancy 



which will be discussed in subsection III B 



The atomic relaxation of the gallium surface antisite 
GaAs is characterized by a downward and inward move- 
ment of the defect with respect to the surface arsenic 
atoms. The equilibrium position is shifted by 0.5 A and 
-0.3 A from the position of the surface anions in the [001] 
and [110] direction. The bond length to the neighboring 



Ga surface atoms is shortened by 3 % compared to the 
defect-free surface. 

Surface antisites on both sublattices do not show any 
charge-transfer levels inside the band gap. The calcula- 
tion of the formation energies of the defects reveals that 
the neutral charge state is always preferable energeti- 
cally. This is in contrast to antisites in the bulk which we 
calculate as double acceptors (GaAs) and double donors 
(AsQa) in agreement with previous theoretical studies 
PI |l6| ]. We therefore predict that surface antisites will 
remain charge neutral on intrinsic as well as on p- or 
n-type doped material and that they do not alter the 
electronic properties of the (110) surface. We note that 
the antisites are electrically inactive at the surface only. 
Once incorporated into the bulk (e.g. at low tempera- 
tures where the surface concentration of defects will be 
frozen during growth) they become electrically active. 
This change in character is rather abrupt as we find anti- 
sites already in the second layer to be similar to bulk ones 
in terms of formation energy and electronic structure in 
agreement with other calculations JIt]]. 

B. Vacancies 

As illustrated in Fig. [I] a vacancy might be derived from 
an antisite by removing the defect atom. By doing so 
the atoms surrounding the defect change their equilib- 
rium positions significantly, and we find a very distinct 
relaxation pattern. The relaxations of the neutral surface 
vacancies on GaAs are shown in Fig.|l]c,d. On both sub- 
lattices they are characterized by a comparatively large 
relaxation of the surface nearest neighbor atoms into the 
surface by 0.5 A (V Ga ) and 0.3A (V As )- The neighboring 
atom in the second layer also relaxes into the void of the 
vacancy by 0.5 A for both defects. The Cih symmetry of 
the defect-free surface is conserved by the relaxation of 
the neutral vacancies. The anion vacancy in the positive 
charge state, however, prefers a nonsymmetric, rebonded 
configuration where one of the surface cations approaches 
the cation in the second layer. The energy gain of this 
configuration compared to the symmetric relaxation is 
O.lOeV on GaP, 0.07 eV on InP J§, and 0.17eV on GaAs. 
The latter value as well as the ground-state geometry is 
in good agreement with the result of 0.16 eV by Zhang 
and Zunger [jl8| for Vas on GaAs. 

Based on our calculations we find that surface vacan- 
cies are electrically active. They have an amphoteric 
character on all three materials investigated, i.e., they 
are positively charged on p- and negatively charged on 
n-type material. Vca on GaAs has two charge-transfer 
levels -&tl(+,0) and £>tl(0, — ) in the lower half of the 
band gap while Va s is a negative U center with -Etl(+, — ) 
in the lower half of the band gap as shown in Fig. |^. We 
can therefore conclude that vacancies are efficient com- 
pensation centers for p- and n-type conditions. 
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C. Adatoms 



As a last type of surface point defects we have cal- 
culated anion adatoms. As illustrated in Fig. |l| we find 
two possible configurations: one in which the adatom 
is bonded to the surface anions (Asn, Fig. Ha) and one 
with bonds to the surface cations (As^, Fig. l]b). In the 
configuration Asn the distance between the adatom and 
the arsenic surface atoms is 2.61 A in the neutral charge 
state while for Asi2 the distance to the surface gallium 
atoms is 2.64 A. The adatom is located 1.1 A (Asn) and 
0.8 A (Ask) above the surface. The defect Asn bound 
to surface arsenic has an electronic structure similar to 
the surface vacancies showing two charge-transfer levels 
£tl(+,0) and £tl(0,-) (Fig.§). In contrast, Asi 2 , 
which bonds to the surface cations, is predicted to be 
negatively charged for all positions of the Fermi level 
except for extreme p-type conditions with -EFcrmi below 
0.1 eV, where the neutral charge state is found to be sta- 
ble. 

Unlike previous first-principles calculations by Yi et al. 
Q the adatom Asi2 is determined as the preferred ad- 
sorption site over a wide range of the surface Fermi level 
(Fig. ||) . Only for p-type conditions near the valence- 
band maximum Asn in the singly positive charge state 
has a lower formation energy than Asi2. For the neutral 
charge states the difference in formation energy is 0.2 eV 
while Yi et al. [[l9| report a difference of — O.leV. This 
difference might be due to the limited size of the sur- 
face unit cell (2 x 2) and the use of the T point only for 
Brillouin zone integration as employed in Ref. E^] . 



IV. FORMATION ENERGIES AND DEFECT 
CONCENTRATIONS 

Based on the formation energies as shown in Fig. [2] we 
identify the dominant point defects on the surface and 
how the concentration depends on the anion chemical 
potential and doping conditions. Under gallium-rich con- 
ditions the gallium antisite GaAs is lowest in formation 
energy (Fig. ||a). It is, however, not electrically active 
and thus can be excluded as a possible dopant or com- 
pensation center. The dominant point defect on GaAs 
under gallium-rich conditions is the arsenic vacancy Va s , 
positively charged under p-type and negatively charged 
under n-type conditions. This conclusion changes if we 
consider arsenic-rich and p-type conditions. Then Vq& 
and Asn, both in the singly positive charge state, are 
almost degenerate in formation energy (Fig.||b). With 
^Fermi in the upper half of the band gap the negatively 
charged As^ is the defect with the lowest formation en- 
ergy among all electrically active point defects. Fig.^Jb 
also shows that the neutral Asc a has a negative forma- 
tion energy under arsenic-rich conditions. This indicates 
that the surface is unstable against the formation of a 
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FIG. 2. Formation energies of surface point defects on 
GaAs(llO) as a function of the surface Fermi level -Bpcrmi- 
(a) gallium-rich conditions with fiAs = — 0.7eV as compared 
to bulk arsenic, (b) arsenic-rich conditions with /xab = OeV. 
Changing to arsenic-rich conditions vacancies on the cation 
and anion sublattice (solid lines) change formation energy 
by —0.7 and +0.7 eV, respectively. Antisites (dotted lines) 
change by +1.4 (GaAs) and — 1.4eV (AsGa) while the adatoms 
(dashed lines) are more preferable by — 0.7eV. 
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complete layer of arsenic. Indeed, first-principles calcu- 
lations of the surface energy show this structure to be 
lower in energy than the cleavage surface assuming ex- 
treme arsenic- rich conditions pp| ] . 

Based on the calculated formation energies we can 
immediately obtain the equilibrium concentration C of 
point defects. Neglecting vibrational entropy contribu- 
tions this quantity is determined by the number of possi- 
ble sites Co where a defect can be formed, the formation 
energy E[ , and temperature T: 

E |(fA 8 - E Fcrmi) 

C = C e . (2) 

Here is Boltzmann's constant. Using this formalism, 
assuming room temperature, and taking into account 
the electrically active defect with the lowest Ef (0.47 eV 
for Vas under gallium-rich conditions with i?Fcrmi at the 
conduction-band minimum) we find a maximum equilib- 
rium concentration of about 10 6 cmT 2 . The concentra- 
tion decreases if we assume a situation where the Fermi 
level is fully pinned at the charge-transfer level of the va- 
cancy at 0.2 eV above the valence-band maximum. Then 
the formation energy becomes 1.27 eV and the concentra- 
tion is negligible (10~ 9 cm" 2 ). A simple analysis shows 
that pinning of the surface Fermi level by defects requires 
rather large concentrations of about 10 12 cmT 2 WA - 
This conclusion is consistent with combined STM and 
photoelectron spectroscopy measurements on InP(llO), 
which showed significant surface band bending at con- 
centrations of phosphorus surface vacancies higher than 
2 x 10 12 cm -2 Q. However, also for the phosphorus sur- 
face vacancy we find an equilibrium concentration (at 
430 K where the experiment has been performed) signifi- 
cantly below the experimentally observed defect concen- 
tration. We therefore conclude that the formation of va- 
cancies is driven by kinetic mechanisms rather than as 
a result of thermal equilibrium. A possible kinetic pro- 
cesses on the surface is diffusion from the bulk driven by 
the lower formation energy of surface vacancies compared 
to their bulk counterparts. 

V. SUMMARY AND CONCLUSIONS 

Density-functional theory has been used to study point 
defects on the (110) surface of GaAs. We identify anti- 
sites on both sublattices to be the defects with the lowest 
formation energy. They are not electrically active and 
can therefore be excluded as compensation centers on this 
surface. Among the electrically active point defects Vas is 
lowest in formation energy under gallium- rich conditions. 
For arsenic-rich conditions and n-type doping Vq & and 
arsenic adatoms are predicted to be most important. The 
calculated formation energies allow the estimation of the 
equilibrium concentrations of surface defects. We have 
shown that the formation energy of all electrically active 



defects is too high to allow for equilibrium concentrations 
large enough to induce significant band bending at the 
surface. The high concentration of defects reported in 
recent experiments is thus related to kinetic effects. 
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